Parity Effect in the Resonant Tunneling 
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A mechanism of the parity effect in the thermally assisted 
resonant tunneling is proposed in the view point of nonadia- 
batic transitions of thermally excited states. In this mecha- 
nism, alternating enhancement of the relaxation is naturally 
understood as a general property of quantum relaxation of 
uniaxial magnets at finite temperatures where appreciable 
populations are pumped up to excited states. It is also found 
that the enhanced sequence depends on the sweeping rate of 
the field. 

PACS number: 75.40.Gb,76.20.+q 



As to the relaxation of metastable magnetization of 
uniaxial nanoscale molecular magnets such as Mni2 and 
Fes, the resonant tunneling phenomena have been paid 
attention and various interesting properties of the phe- 
nomena have been reported flj|-fL0|. The key mechanism 
of the relaxation comes from their discrete energy struc- 
ture due to a finite number of degrees of freedom. The 
eigenvalues of the system are functions of the parameters 
of the system, such as an external field. If we change a 
parameter infinitesimally small, then the system changes 
adiabatically, i.e., if the system is initially in the ground 
state, then it stays in the ground state of the system with 
the current value of the parameter. On the other hand, 
if the changing rate is finite, the system cannot com- 
pletely follow the change of parameter and then so-called 
nonadiabatic transition occurs. For example in uniaxial 
magnetic systems, if we sweep the field very slowly from 
parallel to antiparallel to the initial magnetization, the 
magnetization adiabatically follows the field and reverses 
its direction. This change of magnetization corresponds 
to the tunneling (the adiabatic transition). If the sweep- 
ing rate is fast, then the magnetization only partially 
changes (the nonadiabatic transition). 

In uniaxial magnets, quantum fluctuation is relevant 
only at avoided level crossing points and changes of mag- 
netization only occur at those points. The sweeping rate 
dependence of the probability of staying in the original 
state in this type of nonadiabatic transition has been 
given by Landau [jll), Zener [jl2|, and Stiikclbcrg |ll| 
(LZS). We have studied changes of magnetization in a 
sweeping field from the view point of the LZS mechanism 
|p^-|l6| and proposed to obtain the tunneling gap from 
the magnetization change in a sweeping field JT5[ and also 
explained the step- like magnetization process as a charac- 
teristic feature of nanoscale magnets . The LZS mech- 
anism is pure quantum mechanical and it is independent 
of the temperature. However in experiments, strong tem- 



perature dependences have been observed, which brought 
an idea "thermally assisted resonant tunneling" In 
order to explain this temperature dependence, various 
theoretical attempts have been done Jl7|-|22| . 

The sweeping rate dependences of magnetization pro- 
cess have been also observed in experiments at very low 
temperatures in Mni2 || and also in Fes ||[l0|]- There 
data do not depend on the temperature any more. We 
have pointed out that even in such cases there is still in- 
evitable effect of the environment ]23) , and that we could 
nevertheless estimate the pure quantum transition prob- 
ability. 

Recently new aspects of the resonant tunneling have 
been reported |Q, e.g., the parity effect of the resonant 
tunneling where amount of relaxation changes at the res- 
onant points alternately, and Vi-dependence of the initial 
relaxation of the magnetization at the resonant points 
p4j . In the present Letter, we would like to propose a 
mechanism of the parity effect as a universal property of 
the thermally assisted resonant tunneling. 

In order to investigate characteristics of temperature 
dependence of the resonant tunneling, we have con- 
structed an equation of motion of the density matrix 
where effects of a thermal bath are taken into account 
(the quantum master equation pa]): 



dp(t) 
Bt 



= -i [H, p(t)] - A ([A, Rp{t)) + [A, Rp(t)}^ , (1) 



where TL is the Hamiltonian, p(t) is the density matrix of 
the system and A is a system operator through which the 
system and the bath couple with the constant A. Here we 
set h to be unity. The first term of the right-hand side 
describes the pure quantum dynamics of the system while 
the second term represents effects of environments at a 
temperature T(= /3 _1 ). There R is defined as follows: 



(k\R\m) = ({E k - E m )n {E k - E m )(k\X\m) 
C(w) = — I(— w), and np(u)) = (e 



f) 



where \k) and \m) represent the eigenstates of H with the 
eigenenergies E k and E m , respectively. Here we adopt 
a thermal bath which consists of an infinite number of 
bosons Hb = Y^u^u^' where b u and b]^ are the anni- 
hilation and creation boson operators of the frequency ui. 
We adopt the spectral density of the boson bath I(w) in 
the form I(u>) — Iquj 2 , which is associated with phonon 
reservoir. As to the interaction between the system and 
the bath we adopt a form A ^2 u (b u +&!>)• 

As a more relevant source of a noise at very low tem- 
peratures, we may consider the dipole field from other 
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molecules or/and the hyperfine interaction from nuclear 
spins pfifl . For such noises we have to take into account 
another contribution to R. In this sense the bath treating 
here does not represent the experimental situation very 
appropriately. In the present Letter, however, we discuss 
only general natures which do not depend on the detail 
of the thermal bath. 

Equation ([[]) with X = S x + S z was used to study the 
aforementioned inevitable effects of environments in the 
magnetization process under a sweeping field at a very 
low temperature [ p3[ . For the process the existence of 
interaction is essential but detailed nature of the mecha- 
nism of the dissipation is not important. Thus we could 
discuss the property of the process as a universal prop- 
erty of the relaxation on nanoscale magnets. 

When the temperature goes up and the effect of the 
bath increases, the dissipative process becomes to depend 
on specific features of the bath and the coupling. Thus 
it becomes difficult to treat relaxation without specifying 
nature of the bath. However, in the present Letter, we 
will point out that the aforementioned parity effect is 
a universal property of the resonant tunneling at finite 
temperatures, which is independent of detailed nature of 
the bath. 

Let us consider a general model of an S = 10 uniaxial 
magnet in an external field: 

H = -DSl - TS X - H{t)S z + Q, (2) 

where Q represents extra terms such as (S + ) 4 + (S 1- ) 4 , 
etc. We will propose a mechanism of the parity effect 
as a general property of uniaxial magnets. Thus we put 
Q = 0. 

In Fig. 1(a) we show a magnetization process M(t) = 
(S z )/S for the case with T = 1.0, A = 0.00005, 
T = —0.45. We sweep the magnetic field H(t) — 
ct - H with c = 0.0001 and H = 0.3. In 
Fig. 1(b), the derivative dM/dH{= c~ x dMjdt) is also 
shown. In these figures, we find an alternate change 
of the amount of changes of the magnetization clearly. 
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FIG. 1. (a) Magnetization process and (b) dM/dH for 
c = 0.0001. 
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FIG. 2. Energy structure as a function of the external field 
of the model (||) with D = 0.1, F = 0.45. 

In Fig. 2, energy levels as a function of the field are 
shown. Here we see straight lines along which the mag- 
netization is approximately given by to or to'. These 
lines denoting energy levels for the diagonal parts of the 
Hamiltonian are called diabatic states. At each cross- 
ing point, a small energy gap is created by the off- 
diagonal terms and so-called avoided level crossing is 
formed, where large enhancement of relaxation occurs 
(resonant tunneling). The energy gaps of avoided level 
crossing points are listed in Table I, where we denote 
avoided level crossing points of levels of to and to' by 
(to, to'). There we find that the energy gaps at the same 
horizontal level in Fig. 2 (denoted by the same symbols) 
are about the same. This is easily understood from the 
fact that the gap is of order r'" l ~ m ' where to and to' are 
the magnetizations of the crossing levels for T = [^7| . 

In order to see what processes are going, we show, in 
Fig. 3, the time evolution of distribution of occupation 
probabilities at the zth level, which is expressed by (i\p\i)- 
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FIG. 3. Time evolution of the population (i\p\i). 

This figure shows that the population along the line of 
to = —8 decays at (—8, 5) in Fig. 2 and the populations 
along the line m = —9 and —10 decay at (—9,4) and 
(—10, 3), respectively. These points are shown by squares 
in Fig. 2. 

If the system has an equilibrium distribution in the 
initial state, the population distributes on levels. The 
population at each line decays at an avoided level crossing 
point where the LZS transition probability 

tt(AE) 2 



p = 1 — exp 



2c|? 



(3) 



has an appreciable value. Here AE is the energy gap at 
(to, to'). In Table I the transition probabilities are listed. 
For example the transition probability for c = 0.0001 at 
(to, to') = (—8, 5) (see Fig. 3) is 0.913, while at the point 
(—8, 6) it is 0.07. Thus most of the population of the line 
of to = —8 decays at (—8, 5). The population of the line 
of m — —9 decays very little until (—9, 4) because the 
transition probabilities at (—9, 6) and (—9, 5) are very 
small, i.e., 0.0006 and 0.04, respectively. The population 
of the line of to = —10 decays at the point (—10, 3). 

As we found here, the parity effect simply comes from 
the structure of the energy levels. The structure of energy 
levels in Fig. 2 is inherent to the systems of the uniaxial 
Hamiltonian (0) regardless of the form Q and we expect 
that the parity effect is observed generally in uniaxial 
magnets. 

Here it should be noted that the relevant sequence of 
decays (in the above case H = 0.3,0.5, and 0.7) depends 
on the sweeping rate c. If c decreases, then the transi- 
tion probabilities increase. Thus the populations on the 
lines decay at the circled avoided level crossing points 
before reaching the points of squares. Thus in this case 
relaxation is enhanced at H — 0.2, 0.4, and 0.6 instead of 
H = 0.3, 0.5, and 0.7 in the case of c = 0.0001. In Fig. 4, 
we plot the time evolution of the magnetization process 
in case of 20 times slower sweeping rate c = 0.000005. 
In Fig. 5, we show the time evolution of probabilities at 
levels, where we actually find large decreases of the pop- 
ulation at H = 0.2, 0.4, and 0.6. Furthermore, in case 



of much slower sweeping rate c = 0.00000005, we expect 
that the transitions occur at the points of triangles in Fig. 
2 (see also Table I). We do not demonstrate it because 
it takes 100 times longer simulation time. This sweep- 
ing rate dependence is also a general property of uniaxial 
magnets and we expect the shift of the sequence to be 
also found in experiments. 
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FIG. 4. (a) Magnetization process and (b) dM/dH for 
c = 0.000005. 
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FIG. 5. Time evolution of the population (i\p\i). 

Finally we would like to point out a strange property 
of the system when the Hamiltonian (|J) includes the in- 
teraction 



Q = C[(S + r + (S-) 



(4) 



which has been discussed in literatures As has 

been pointed out by Wernsdorfer and Sessoli M, in the 
present model energy gaps at the avoided level crossing 
point change nonmonotonically with the value of C {V 
is fixed) and even gapless points exist, which causes ir- 
regular behavior of resonant tunneling and simple parity 
effect is disturbed. The reason why gapless points appear 
is not clear at this moment, which would be an interest- 
ing problem. 
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from the Ministry of Education. 
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H 


AE 


(m, m') 


p(c = 0.00000005) 


p(c = 0.000005) 


p(c = 0.0001) 


0.000 


0.1554D-05 


(-10,10) 


0.3795D-05 


0.3795D-07 


0.1897D-08 


0.000 


0.1402D-05 


(-9,9) 


0.3432D-05 


0.3432D-07 


0.1716D-08 


0.000 


0.8514D-04 


(-8,8) 


0.1413D-01 


0.1423D-03 


0.7116D-05 


0.000 


0.9031D-02 


(-7,7) 


0.1000D+01 


0.8396D+00 


0.8745D-01 


0.100 


0.1466D-05 


(-10,9) 


0.3555D-05 


0.3555D-07 


0.1778D-08 


0.100 


0.5161D-05 


(-9,8) 


0.4923D-04 


0.4923D-06 


0.2461D-07 


0.100 


0.9741D-03 


(-8,7) 


0.8629D+00 


0.1968D-01 


0.9932D-03 


0.101 


0.5483D-01 


(-7,6) 


0.1000D+01 


0.1000D+01 


0.9735D+00 


0.200 


0.1391D-05 


(-10,8) 


0.3375D-05 


0.3375D-07 


0.1688D-08 


0.200 


0.7220D-04 


(-9,7) 


0.1018D-01 


0.1023D-03 


0.5118D-05 


0.200 


0.8156D-02 


(-8,6) 


0.1000D+01 


0.7752D+00 


0.7191D-01 


0.300 


0.3558D-05 


(-10,7) 


0.2340D-04 


0.2340D-06 


0.1170D-07 


0.300 


0.7418D-03 


(-9,6) 


0.6841D+00 


0.1146D-01 


0.5760D-03 


0.302 


0.4497D-01 


(-8,5) 


0.1000D+01 


0.1000D+01 


0.9131D+00 


n /inn 






U.ooZOJJ — uz 


n 33^1 "n c\a 

U.OOOl-L/ — U4 


n i RfifiTi n^ 


0.400 


0.5846D-02 


(-9,5) 


0.1000D+01 


0.5356D+00 


0.3762D-01 


0.500 


0.3846D-03 


(-10,5) 


0.2664D+00 


0.3093D-02 


0.1549D-03 


0.502 


0.2909D-01 


(-9,4) 


0.1000D+01 


0.1000D+01 


0.6404D+00 


0.600 


0.2938D-02 


(-10,4) 


0.1000D+01 


0.1761D+00 


0.9639D-02 


0.702 


0.6935D-01 


(-10,3) 


0.1000D+01 


0.1000D+01 


0.9970D+00 



TABLE I. Energy gaps at avoided level crossings and LZS probabilities 
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